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Near-threshold production of charged boson pairs is considered within the framework 
of the model of unstable particles with smeared mass. The results of calculations are 
in good agreement with LEP II data and Monte-Carlo simulations. Suggested approach 
significantly simplifies calculations with respect to standard perturbative one. 
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1. Introduction 

The measurements of VT-pair production at LEP II provided us with an impor- 
tant information about the mass of W boson and non-abehan triple gauge-boson 
couplings. To extract the exact information from VF-pair production we have to 
calculate the radiative corrections (RC's), which give a noticeable contribution 
to the cross-section. Ideally, one would like to have the full RC's to the process 
e^e~ — > W^W~ — !■ 4/. In practice, this problem is very complicated and can not 
be considered analytically. For discussion of the LEP 11 situation and strategy it is 
useful to distinguish three levels of sophistication in the description of the M^-pair 
production [HE]: 

1) On-shell VF-pair production, e+e~ W^W~ , with consequent on-shell W 
decays. All 0{a) RC's to these processes are known. 

2) Off-shell production of W pairs, which then decay into four fermions. Full set 
of RC's is very bulky for the analytical observation and analysis. 

3) Full process e+e^ 4/ with an account of the complete set of the 0{a) 
corrections. This problem leads to the additional diagrams with the same final 
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states, and complete electroweak 0{a) corrections are described by many thousands 
diagrams. 

On-shell W^-pair production was considered in Refs. [Hill [3], where the cross- 
section of the process e+e" W^W~ was given. At tree level, this process is 
described by three diagrams shown in Fig. 1. The complete 0{a) radiative correc- 
tions, comprising the virtual one-loop corrections and real-photon bremsstrahlung, 
were calculated and represented in Refs. [4]- [11]. The description of the on-shell 
VF-pair production and consequent decays with an account of RC's was fulfilled in 
Refs. [12]-[IH]- Off-shell production of VF-pairs, which then decay into four fermions, 
was considered in Ref. |19,. 

In description of the W- and Z-pairs production we should take into con- 
sideration the fact that the gauge bosons are not stable particles and the real 
process is not e+e~ W^W^ , ZZ [2\. This is only an approximation with a 
level of goodness, which may depend on several factors, while the real process is 
e+e" W~^W~ , ZZ 4/. There are many papers, devoted to comprehensive 
analysis and description of all possible processes with the four-fermion final states. 
Because of a large number of diagrams, describing these processes, the classifica- 
tion scheme was applied in Refs. [25-[23]. The possible processes are divided into 
three classes: charge current (CC), neutral current (NC) and mixed current (MIX). 
Born processes e+e~ W^W~ ^ZZ are designated as CC03 and NC02, which 
correspond to three charge current and two neutral current diagrams. According to 
this classification the off-shell VF-pair production with consequent W decay can be 
described in the framework of the Double-Pole Approximation (DPA) [53]- [55]. The 
DPA selects only diagrams with two nearly resonant W bosons and the number of 
graphs is considerably reduced [23] . 

Complete description of the total set of 4/-production processes with an ac- 
count of RC's is not analytically available due to a huge number of diagrams 
and presence of non- factorable corrections. But the complete EW 0(a) corrections 
have been calculated for some exclusive processes, for instance, for the processes 
e+e" VtT^ [i~ v^j^^ udyTVp,^ and udsc [27l|28]. Because of complexity of the prob- 
lem, some approximation schemes are practically applied, namely. Semi- Analytical 
Approximation (SAA) [2l I29j. improved Born approximation [30] . an asymptotic 
expansion in powers of the coupling constant of the cross-section [31| , fermion-loop 
scheme, etc. (see Introduction in Ref. [27] [28]). There are many computer tools 
of calculations, for instance, Monte-Carlo (MC) simulations, such as RacoonWW 
[28l|3l|33] and YFSWW [Ml [35l|36] . All above mentioned methods are based on 
the traditional quantum field theory of unstable particles (UP's) [2]. At the same 
time, there are some alternative approaches for description of the UP's, such as the 
effective theory of UP's [33-[39] and the model of UP with smeared mass [40 ] I42 j . 

In this paper, we suggest the description of the near-threshold ly-pairs produc- 
tion within the framework of the model of UP with smeared mass initially pro- 
posed in Ref. [40]. The model is based on the time-energy uncertainty relation 
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AE • At ~ 1 (c = ?i = 1). It follows from the equation of motion in the Heisen- 
berg representation which describes the evolution of the non-stationary quantum 
system [43]. In the case of the unstable particles, At is the lifetime and AE is the 
value of the mass smearing Am in the rest-frame system |40l |43l [44] . In the model 
under consideration the UP is described by a state with smeared (fuzzed) mass in 
accordance with the uncertainty relation. So, the processes e+e~ W~^W^ , ZZ 
are described in a traditional way, i. e. in a stable particle approximation, but the 
phase space is calculated for the states on the smeared mass-shell. In the frame- 
work of the model, full process e+e~ W^W~, ZZ 4/ is divided into two 
stages e+e~ W^W'^ , ZZ and W^W~ , ZZ 4/ due to exact factorization 
at tree level [45]- [47]. For description of the bosons in the final state we use the 
model polarization matrix which differs from the standard one |42[ 147] (see the next 
section). 

The model was applied for description of the Finite- Width Effects (FWE's) in 
many low- and high-energy processes involving the UP with large width [42], |45] - 
[48j . In particular, the approach was successfully applied to the process e+e~ ZZ 
in Ref. [42]. In this paper it was shown that the results of the model are in good 
agreement with the LEP II data and turned out to be very close to the corresponding 
results of MC simulations. So, it is reasonable to apply the same method for the 
case of VF-pair production. 

It should be also noted that the model under consideration directly leads to 
calculation schemes which are in close analogy with such standard approaches as 
the Convolution Method (CM), Narrow- Width Approximation (NWA) and SAA. 
However, the model treatment has some noticeable distinctions which are discussed 
in detail in Refs. [IS]- [33 (see also the next section). The principal distinction 
between the standard and model treatment of the FWE's takes place in description 
of UP with large width [57] and the mass splitting in the neutral meson systems 

m- 

The paper is organized as follows. In the second section, we give a short de- 
scription of the model and define the status of our calculations. The cross-section of 
the process e+e^ W^W^ at tree level is derived in the framework of the model 
[40] in this section. Section 3 contains the calculation strategy with taking into ac- 
count of the radiative corrections. In this section we also represent the results of our 
calculations, MC simulations and LEP II data. Some conclusions concerning the ap- 
plicability of the method were made in the last section. We note that the aim of this 
investigation is to test the model approach for the case of near-threshold boson-pair 
production and analyze the possibility of its improvement by an accounting of the 
radiative corrections which are not included into the effective field of UP. 
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2. The model cross-section of the near-threshold M^-pair 
production at tree level 

Firstly, we give a short description of the model of UP with smeared mass. The 
model wave function of the UP is 

/ <^a{x,^l)uJ{^x)dn, (1) 



where $a(a^, is the standard spectral component, which defines a particle with a 
fixed mass squared = in the Stable Particle Approximation (SPA). The weight 
function lo{^) is formed by the self-energy interactions of UP with vacuum fluctua- 
tions and decay products. This function describes the smeared (fuzzed) mass-shell 
of UP. Thus, the smearing of mass is caused, on the one hand, by instability accord- 
ing to formal uncertainty relation and, on the other hand, by stochastic interaction 
of UP with the electro- week vacuum fluctuations [46] . 

The commutative relations for the model operators have an additional (5-function 

fc(fc,/i), mq,^Ji')]± = 5{^i-^i')6{k-q)6^p, (2) 

where subscripts "±" correspond to fermion and boson fields, respectively. The 
presence of S{fi — n') in Eq. ([2]) means the following assumption: the acts of creation 
and annihilation of the particles with various /i (the random mass squared) do not 
interfere. Thus, the parameter fi has the status of physically distinguishable value 
of a random rri^. 

The model Green functions for the vector and spinor fields in momentum rep- 
resentation have the convolution form: 

'-^ — — p{^^)<^^^ , (3) 

- 11 + le 

and 

m-^I^^^^Pi,)d,, (4) 

where p{fi) = |w(^)p. 

Further, we consider the model amplitude for the simplest processes with UP in 
the initial or final state. The expression for scalar field is 

= J ^ip)dt^ J ^^^^e^^'^^dq, (5) 



where gJJ — y^cp + n and a^(g, /i) are the creation or annihilation operators of UP 
with the momentum q and mass squared = /i. Taking into account Eq. one 
can get 

[a-(fc,M),0+(:r)]-, r(a:),a+(fc,/i)]-- ±a.. ^ (g) 

(27r)3/2 - - 



19, 2009 6:15 WSPC/INSTRUCTION FILE WWcrsecl 



Near-threshold W-pair production in the model of unstable particles with smeared mass 5 



where = ^/k^^+Jt. The expressions ^ differ from the standard ones by the factor 
uj{n) only. From this result it follows that, if d"'"(fc, ^)|0) and {0\a~ {k, fj,) define UP 
with the mass m — y//I and momentum k in the initial or final states, then the 
amplitude for the transition $ is 

A{k,^i)^u;{^l)A''{k,^i), (7) 

where A''*(fc,/i) is the amplitude in the SPA. This amplitude is calculated in the 
standard way and can include the higher corrections. Moreover, it can be an effective 
amplitude for the processes with hadron participation. From Eq. ([7]) it follows that 
the differential (over /x) probability of transition is dP{k,fi) = p{fi)\A{k, fi)\'^d^. 

To define the transition probability of the process $ —> (jxpi , where (j> is UP with 
a large width, we should take into account the status of parameter /i as a physically 
distinguishable value, which follows from Eq. ([2|). Thus, the differential (over A:) 
probability is 

dT{k) ^ J dr'*{k,fi)p{n)dfi. (8) 

In Eq. ([5]) the differential probability dr^'^{k, /i) is defined in the standard way (the 
SPA). 

If there are two UP's with large widths in the final state of decay <f> 0102, 
then in analogy with the previous case one can get the double convolution formula: 



r(m$) = J y r^*(m$;/xi,/X2)pi(Aii)p2(Ai2)dMi<^Ai2 • (9) 
The polarization matrix for the case of vector UP in the final state has the form 

e 

In the case of spinor UP in the final state we have 



1 



(9)V (9)"^(^TVM)a/3, (11) 

1/ 



where the summation over polarizations is implied and = y^q'^ + ii. The same 
relations take place for the initial states, however, one have to average over their 
polarizations. 

The most important element of the model is the probability density p{fi) which 
describes the smearing of UP mass. The various definitions of p(/x) were discussed 
in Ref. [46], where the Lorentzian (Breit-Wigner type), Gaussian and phenomeno- 
logical distributions have been considered. The Lorentzian distribution was derived 
by matching the model propagators ([3]) , ([4]) and standard dressed ones in the Breit- 
Wigner form. 

One of the important properties of the model is the exact factorization of the 
processes with the UP in an intermediate state. In the frame of the effective theory 
of UP, which follows from the model, the factorization leads to the convolution 
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formula for the decay rate [ITl and factorized formula for the cross-section |45) . 
These results are derived by straightforward calculations at tree level without any 
approximations. So, the model provides the formal basis for the CM, SAA and 
NWA, which are the approximate approaches in the traditional treatment. The 
generalization of the factorization method to the complicated processes of scattering 
and decays with two or more UP in intermediate states was considered in Ref. [47] . 
This method can be applied for description of the boson-pair production and decays 
in the factorized form for double-pole set of diagrams (see comment to tree-level 
result at the end of this section). 

In this work, we use two principal elements of the model [40|,l46j - the convolution 
structure of the transition probability (an analog of the expression jSl)) and the 
polarization matrix for UP in the final state pO|l . Using these expressions, we get 
the model Born cross-section of M^-pair production in the following form 



'WW 



f^wwi^'^ /^i,M2)/Oi(A^i)/02(Ai2) d^ii , 



(12) 



where a^^{s; Hi, fJ,2) is the Born cross-section which is calculated in the standard 
way for fixed bosons masses /^i = ml and /i2 — ^2 (SPA). 




V >r 



w 



+ 




Fig. 1. Feynman diagrams for the process e+e — > W^W 



The Born cross-section is defined by the sum of two diagrams shown in Fig. 1 
and can be represented as 



7ra 



128s sm Ow 



(13) 
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where dimensionless function F{s; xi, X2) is defined by the expression 
16 

F{s; xi,X2) = , r^{3(a2 _ ^2)(^2 _ ^2 ^ 2(1 + a))(l - xzfL{a, b) 

6(a^ — b^)(l — xz) 

+ a) + a (4 + a))(l — xz)L{a, b) 
+ 2A(a, b){2b^ - 3a^ - 10a - l){b^{l - 2xz) - a(l - 3xz) ~ xz)\ 
+ A(a, h)\x\}?{a, b) cos(46'vi/)(26^ - - 10a - 1) + \2c?z\ 
- a^{Zb^(3x\ - 2xz + 1) - 49x| + 30xz - 15) - 2a{b^{\'dx\ - lOa^z + 5) 
+ %x\) + 26'*(3x| - 2xz + 1) - 2b^{lx\ - ^'^xz + 8) - 2.t|]}. (14) 

In Eq. (|14p the dimensionless variables a, 6, xi, a;2, a;z and the functions L{a,b) and 
A (a, 6) are defined as follows 

1 — a — A(a, 6) 



i(a, b) = In 



1 — a + A(a, b) 



X{a,b) = \J\ - 2a + 62, 



xi 9 = , a = a;i+a;2, = Xi-X2, xz = ■ (15) 

s s 

With the help of the expressions (fT3|) - p^ the model Born cross-section is repre- 
sented in the following convolution form 

^ww{s)^— — . 4^ / dxip{xi,s) p{x2,s)F{s;xi,X2)dx2, (16) 

128s sm Ow Jo Jo 

where in analogy with the case of Z-pair production [42j we use the redefined 
dimensionless probability density of Lorentzian type 

1 G{x,.3) V^rt?-* 3ax 

P(X, S) = ^2 I n2/ T' ^) = = , . 2a ' i^^) 

IT (x — xw) + Cjf (x,s) s 4sm^i^ 

where x = p/s, xw = M^/s. The expression (I16p turns into the standard expres- 
sion for the on-shell cross-section cr^^y(s), when p(/i) — > (5(/i — M^), i.e. in the 
limit of fixed W mass. 

In Fig. 2 we represent the Born cross-section in the fixed-mass approach (solid 
line) and in the smeared-mass approach (dashed line). Besides, in this figure we 
show the experimental LEP2 data in order to illustrate the necessity of radiative 
corrections. One can see that the model approach leads to the smearing of the 
threshold, that is to the result which is similar to the standard one with an account 
of the Finite- Width Effects (FWE). In other words, the model description of boson- 
pair production is in close analogy with the standard description of the off-shell 
boson-pair production. Moreover, the model convolution representation of the cross- 
section in the form ([T2|) or ([T6|) is formally similar to the SAA [21 [21]. However, the 
SAA is constructed as an approximation for the exclusive process like e^e^ — > 
W^W~ — > /i 72/3 A, which then is generalized to the inclusive process with full 
set of the 4/ final states. This approach is based on the approximate factorization 
of the total cross-section a{e+e- /1/2/3/4) — > cr(e+e- W+W-)Bt{W+ 
/i/2)Br(W^~ f-ifi). In the frame of the model [40l|46], the factorization is exact 
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160 170 180 190 200 210 

4^, GeV 

Fig. 2. On-shell (solid line) and model (dashed line) Born cross-section of the process e+e~ — » 
W+W-. 

(see also Refs. [13 U?]), and the expression (fT^ can be directly derived [37] for 
the inclusive process e^e" — s- 4/ in the double-pole approach without any 
approximations. 

From Fig. 2 it follows that the use of the effective model fields, which describe 
the UP with an account of the self-energy type corrections, is not sufficient. We 
have to take into account the rest radiative corrections for the realistic description 
of the measured cross-section. 

Now we estimate the uncertainties of the model calculations at the effective tree 
level, which are caused mainly by the definition of the function /o(/i). As was shown 
in the framework of the effective theory of UP [45J-[47|, this function results from 
the factorization of full process of production and decay of UP. Let us define the 
uncertainty as a deviation of the model calculation from the standard one. The part 
of the model amplitude which describes the decays W — > Ivi is 

The standard expression for the amplitude follows from Eq. (fT5)) after the change 
— ^ in the numerator of the dressed propagators of unstable bosons W. We 
assume that the denominators P{q^) in both cases are the same (in the Breit-Wigner 
or complex pole form). As was mentioned above, the use of the standard propagators 
does not lead to the exact factorization even at the tree level. In the standard 
approach this effect takes place in the Narrow- Width Approximation (NWA) , while 
in the framework of the model under consideration the factorization is exact due to 
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specific form of propagator's numerator |47j . Using the equality 

lqa{i ~ To)l^a = mj{l - j5)Va, (19) 

where a = 1,2, we get 

|M™°'*|2^1 + 2— + 2^, (20) 

11 92 

Where q = y/(q ■ q). The same expression takes place for the standard amplitude 
squared after the change qa Mw- As a result, we have the relative devia- 

tion of the model partial cross-section from the standard one (mi ~ m2 = mf)'- 

e,^A^[l-M^ r ^dq% (21) 

Mw Jmj q 

From (|21[) with the help of the Breit-Wigner approximation for the function 
piq"^) we find that the maximal deviation is for heavy fermions, for instance, for the 
T-lepton pair and b, c-quark pairs. However, in the last case the e is suppressed by 
small CKM elements |C/cbP- Thus, from ((2T|) it follows that e^ax = Cr ~ 10~^, that 
is the uncertainty of our approach is an order of 0.1%. From this simple analysis 
we make a conclusion that an error, which is caused by the model approach at tree 
level, noticeably less than 1%. So, the main uncertainty can be caused only by the 
implantation of the radiative corrections into our scheme of calculation (see the 
next section). 

3. The model cross-section of W-pair production with radiative 
corrections 

In this section we discuss the strategy of the RC's accounting and represent the final 
results of calculations. As it was shown in Refs. [l5lll6|, the model description of UP 
is equivalent to some effective theory of UP, which includes the self-energy type RC's 
in all orders of perturbation theory. Moreover, the UP is the non-perturbative object 
in the vicinity of the resonance. So, the traditional program of RC's calculation is 
not valid in the framework of the model. We have no well defined set of the diagrams 
which is gauge invariant and renormalized. The model of UP [40l |46] is effective and 
not gauge one, and we have no any rigid criteria for definition of such a set. So, we 
keep the strategy which is based on the simple phenomenology and was successfully 
applied in the case of Z-pair production 42J . 

We do not take into account any corrections to the final states W, because of the 
effective nature of these states in the framework of the model. We use the effective 
couphng a{Mw) = 1/127.9 in the vertex with the final W^-states and a — 1/137 
in the RC's. So, the principal part of the vertex corrections is effectively included 
into the coupling, and the low-energy behavior of the bremsstrahlung and radiative 
corrections to the initial states is taken into consideration. 

The set of corrections, caused by the final state interactions in the two s-channel 
diagrams in Fig. 1, is included into the effective coupling a(Mn/). The principal 
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part of the so-called Coulomb singularity contributions, which were considered in 
Refs. [1], [27] and [3Tj, can be also absorbed by the effective coupling. The one-loop 
calculation shows that this correction gives from 5.7% at the threshold to 1.8% at 
190 GeV [1], while the total change of the effective coupling a{Mw) with respect to 
a is near 7%. In the calculation we explicitly take into account the 0{a) corrections 
including soft and hard bremsstrahlung, which are not described by the model and 
by the effective coupling. The real and virtual electromagnetic radiation should 
enter into the set of these RC's and mutually compensates the total IR divergences. 

The program of RC's calculations, which is similar to above discussed one, was 
fulfilled in the series of papers (see, for example, Ref. [11] and references therein) 
for the case of the on-shell VF-pair production (the limit of fixed masses /ii = /i2 
M^f). The analytical expression for these corrections is represented in compact and 
convenient form in Ref. [11] . We generalized this expression to the case of smeared- 
shell VF-pair production, that is for arbitrary values of mass parameters /ifc, and 
applied it in our calculations. As a result, we get the cross-section (Tvvw(s; /^i, M2) 
for the case of VF(^i) and W{ii2) production including above described corrections 
in the following form (see also Ref. [TT] ) 

(TWw{s■,^il,^JL2) ^ / p^(fc)cr^^y(s(l - fc),^l,^2)rffc, (22) 





where p^{k) is the photon radiation spectrum [49]-[5T], k = E^/Ei, is the photon 
energy in units of beam energy and s(l — k) is the effective s available for the W- 
pair production after the photon has been emitted [TI]. In the case of the on-shell 
VF-pair production (/^i — ^2 — ^^w) ^^'^ value kmax = 1 ^ 4Af^/s is the maximal 
part of photon energy. The generalization of this value to the case ^ fi2 leads to 

krnax = l-2 \ = A (^i , /i2 ; s) . (23) 

s 

The photon distribution function is written in the form [llj 

p^ik) - (3kf'-\l + 6^,+' + ...) + S'l + (24) 

where we keep 0{a) corrections only (i.e. (5„>i — 0). The corresponding corrections 
are given by (u + s = virtual+soft, h ~ hard) [TT] : 

P =—iL-l), L = ln^, 

- -(|i + ^ - 2), S'^ = -(1 - i)(2 - k). (25) 

Finally, we get the corrected expression for the cross-section of PF-pair produc- 
tion in the form 

awwis) ^ .4 / dxi / dx2 / rP7(fc) 

128s sm Ow[Mw) Jo Jo Jo 

p{xi, s(l - k))p{x2, s{l - k))F{s{l ~ k);xi,X2). (26) 
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where the functions F, p, were defined before in Eqs. (dH), (HH), and we also 
take into account the effective QCD correction factor kqcD = 1 + O.ISS/tt [52] . 




T-*-T-* 



160 170 180 190 200 210 

V^, GeV 

Fig. 3. Model (dashed line) and Monte-Carlo RacconWW and YFSWW (solid lines) cross- 
sections of the process e+e~ — > W^W~ . 

The model cross-section aww (s) was calculated numerically and represented in 
Fig. 3 as a function of s by dashed line. The results of MC simulations, RacconWW 
[32l [33] and YFSWW [34l |35] , are represented for comparison by two barely dis- 
tinguishable solid lines, and the experimental LEP II data [53] are given with the 
corresponding error bars. From Fig. 3, one can see that the model cross-section with 
RC's is in good agreement with the experimental data. Moreover, the deviation of 
the model from MC curves is significantly less then the experimental errors (< 1%). 

In the previous section, we have got an estimation of the tree level uncertainty 
which turned out to be significantly less than the value of radiative corrections. 
So, the total uncertainty of the model approach mainly depends on the set of the 
RC's which we take into account. From the above described strategy of the RC's 
accounting, it follows that the principal value of the error can be caused by the 
part of the Coulomb corrections, which follows from the box diagram [1] and gives 
most likely less then 1%, and by the non-factorable corrections, which destroy the 
convolution structure of the total cross-section. A comparison between the DPA and 
the predictions based on the full 0{a) corrections reveals differences in the relative 
corrections < 0.5% f27 and 0.9% f5T. 

From the results of our calculations it follows, that the model approach provides 
the accuracy which is sufficient for the LEP II data description in the near-threshold 
energy range. Besides, the contribution of the non-factorable corrections in the 
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cross section is less than the experimental errors. However, the value of all non- 
considered corrections can be maximally an order of LEP II uncertainty of the total 
cross-section, and this point needs an additional consideration. Moreover, the actual 
status of the calculations concerns rather the testing of our approach than the tool 
for precise investigations. But, we believe that the approach due its simplicity and 
physical transparency can provide the basis for construction of such a tool. 



4. Conclusions 

A large number of the first order Feynman diagrams, which contribute to the pro- 
duction of four fermions in e+e" interactions, depends on the specific final states. 
So, the detailed classification was suggested for the description of these processes. 
The most important request for WW physics concerns the 0{a) radiative correc- 
tions in the DPA. Inclusion of the complete EW corrections significantly complicates 
the calculations which became not available in the analytical form in the case of the 
full set of 4/-processes. So, the various approximation schemes have been worked 
out together with development of the MC simulations. 

In this paper, we applied the model of UP with smeared mass for description 
of the W-pair production. The model describes the process e+e~ W~^W~ as 
H^-pair production, where W's are on the smeared mass-shell. This approach is 
similar to the standard description of the off-shell VF-pair production in SAA. We 
have taken into account the soft and hard initial state radiation and a part of the 
virtual radiative corrections which are relevant in the framework of the model. 

It follows from our results that the model is applicable to description of the 
near-threshold boson-pair production with LEP II accuracy. We get the total cross- 
section which is in good accordance with the experimental data; it coincides with 
the MC calculations with a high precision. At the same time, the model provides 
a compact analytical expression for the cross-section in terms of convolution of the 
Born cross section with probability densities (or mass distributions) of W bosons. 
However, we did not fulfill the detailed analysis of an accounting of the EW cor- 
rections, so this rather phenomenological formalism can not be directly applied for 
the precise description of the boson-pair production at high energies and for future 
experiments at ILC. It is reasonable to consider the possibility of improvement of 
the approach and its applicability at the energies far from the near-threshold range. 
We leave this analysis for a separate study. 
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